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Abstract 

A New Design/Analysis system for the flews 
through tiirborachirkery is currently being developed 
for studying turtxarachinery problems with an 
axisy^aextric viscous/ invisc id "average-passage" 
thrcucdiflcw ccx&. Th& advantage of this approach, 
c gpred to streamline curvature codes, is that the 
solutions obtained simulate sane of the 
unsteadiness, conpressibility and viscous effects of 
a multistage turtxrachine . The Design/Analysis 
system consists of three elemental parts, the 
axisymmetric block grid generator, the blade surface 
element code, and the axisymmetric flow code. Each 
element of the system will be discussed and the flow 
solutions for three axisynmetric geometries will be 
shown compared to experimental data where available. 
The canpsitations are shewn to be in very good 
agreement with test data for SR7 Spinner Body and 
Transonic Boattail Q&xn&txy obtained in the wind 
tunnels at NASA Lewis Research Center. While the 
VIADAC Rotor 67 Fan results were compared to PARC2D 
calculated results and shewn to be in very good 
agreement. 

Introduction 

Early, the design methodology used for 
multistage turbcarachinery has relied heavily on the 
streamline curvature approach [1-7] for confuting 
the flows through modem gas turbine engines. These 
systems have their foundations in the early 
discussions of Wu [1] and Marble [2] on the three- 
dimensional flows in turbcnachinery. Even though 
these methods have proven to be quite successful, 
the current approach has been to model the flews 
with Euler and Navier-Stokes formulations. This 
approach has been taken in the develcpment of a new 
design/analysis system for flews through 
turbciiachinery at NASA Lewis Resea r ch Center. 

Smith [3], Novak [4], and Jemions and Stow 
[6,7] have all described in detail, the advantage of 
using a quasi-three-d Imensicnal streamline curvature 
approximation to design. These formulations are 
based upon Wu's [1] work with S 1 and S 2 
streansurfaces. The thrcughflcw calculations are 
done on the S 2 streamsurface, while the blade- to- 
blade cenputatuons are done an the S, streamsurface. 
In general , the axisynmetric ccxiputaticns are 
calculated with sane built in enpixic&l data base to 
account for secondary flews and spanwise mixing 
effects [8-12] . This methodology has been #£te 
successful in cxxrputing the performance and design 
features of the blades for turbcinachinery . 

More recently, Ni £13] and Holmes [U] have 
described methods of axputing the 3-D inviscid 
flews for isolated blade rows in tnr ba^rhinery, 
and Hah [15] and Chima [16] have developed methods 
of computing 3-D Viscous flows for simgle-hlade-rw 


turbcnechinery prolans. One limitation to all of 
these methods is that they are primarily used for 
analysis of single-blade-rows and not applicable to 
the design of turbonachinery. 

Dawes [17], and Denton [18] have preposed 
various methods of calculating the 3-D flows in 
multi-blade-row turbomachinery by P^ing 
tangentially average flow quantities between blade 
rows in the computation. In Denton's method, the 
flow field effects of the adjacent blade row are 
axisymnetrically averaged at an ^ axial sts. xcs\ 
approximately mid-way between the blade rews. 1, 
flow calculations neglects ary unsteady effects 
which might arise from the adjacent blade rows. The 
approach by Dawes is very similar to that of Denton, 
however, Dawes allows each blade row to be 
calculated either fully 3D, or axisyirretrically with 
foro^, loss and deviation either fran a 




Adairczyk's "average-passage'' approach [19], 
describes a model in which the 3D computations are 
performed on each individual blade row with the 
results fron adjacent blade rows being represented 
as force terms. 3D solutions obtained fron 
individual blade rows are averaged circumferentially 
to cbtain the force terms neede d to perform the 3D 
calculations on the other blade rows. This method 
nofels the deterministic unsteadiness of the flow 
field through a multi-stage turbanachine in a time- 


This "average-passage" approach is being 
used in the development of the new design/analysis 
flow solver. The flow code being developed at NASA 
Research Cen ter is similar to the traditional 
quasi— 3D design codes. The flow solver is 
axisymnetric and can be run inviscidly with assumed 
or calculated blockages, or with the visocus terms 
concerted. The blade forces for each blade row are 
frtxn blade- to-blade solutions, correlated 
data or force ncdel, or fron a full 3D solution, 
aiade-to-blade amputations can be obtained with 
existing or with the blade-to-blade solver 

provided for the design system. This method 
represents an approximation to the full 3D 
amputations. With the proper modeling of the 
spanwise mixing, secondary flow effects, and blade 
force terms, this model should closely represent 3D 
flows through multi-stage taxfeoraachinery without 
actually having to compute the full 30 solutions. 


The Desion/Analvsis Syste m 


The design and analys 
currently under daveloptaent a 
Oenter have been separated 
elements. These areas ccnsist 
interactive grid generator 
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restart code, TTGGERC, the interactive blade element 
geometry generator, IBBGG, and the viscous/irrviscid 
rnulti-blade-row average passage flew solver, VIADAC. 
The block diagram, shewn in figure 1, provides a 
flew chart as to hew information is shared between 
codes. 



Figure 1. Flow Chart for NASA LERC Turbomachinery 
Design/Analysis System 


The large blocked areas represent the 
activities currently being undertaken at Lewis 
Research Center. The other areas, such as the 
axisymmetric streamline curvature solution, quasi-3- 
D blade-to-blade solution and full 3-D solution, 
have already been developed as separate flew 
solvers. A quasi-3D blade-to-blade solver is being 
provided to compute the blade forces if no other 
method of obtaining these forces is available. The 
codes are modular and the design/analysis system can 
be used with existing codes. 

T\rrbcmachiner/ Interactive Grid Generator Energy 
Restart Code (TIGGERC) 

One of the key elements in obtaining the 
axisymmetrie flow solution is the grid generation. 
An interactive block grid generator, called TIGGERC, 
has been developed for any type of txutxriachinery 
problem. The grid generator creates an axisynmetrie 
mult i-bl ock grid for jinyjduct , single rew or multi- 
blade- row txudxxoaciunery problem. 

TTGGERC is fully interactive and can be used 
to modify flcwpath coordinates and grid packing. 
The code can be run on the Silicon Graphics 4D/GT 
Personal Iris or any other Silicon Graphics 
Workstation. The input coordinates can be entered 
interactively or by an input data set. TTGGERC is 
mouse driven and several menu features have been 
added to make the system as general as possible. 
For example, Bezier curves are available to modify 
block boundary coordinates. Grids are generated 
using a hyperbolic tangent or algebraic distribution 
of points on the boundaries and the interior points 
of each' block are distributed using a trans finite 
interpolation. Primarily, TTGGERC can generate a 
Blocked H-Grid for the axisynmetric code like the 
examples of a small multi-stage engine and NASA 
Rotor 67 Fan stage, shown in Figure 2. However, 
The grid generator can also be used to generate C- 
Grids, I -Grids or O-grids, as shown in Figure 3. 


The code is general enough to handle any arbitrary 
geometry for 2D grids. 

TTGGERC can be used to sequentially generate 
several different grids on the same gecretry. This 
feature provides the ability to generate inviscid or 
viscous grids, repack the grid, and add more points 
along the surface, all without restarting the code. 
TTGGERC is to be linked to the Interactive Blade 
Element Geometry Generator by the leading and 
trailing edge locations of the blade, blade 
stacking, number of blades and blockage. A more 
detailed discussion on the menu system used by 
TTGGERC and the internal workings of the grid 
generator will be published later this year. 


interactive ^ade Element Geometry Generator (IBBGGj 

The second key element for the design system 
is the interactive blade element geometry generator, 
IBB GG. This code is a crucial part of any 

design/analysis system. The blade element code 
computes the upper and lower blade surface geometry , 
stacks the blade along a prescribed stacking axis, 
provides choke margin information for cefnpressors , 
and allows the chord, solidity, thickness 
distributions, blade angle distribution or blade 
surface distributions to be altered interactively by 
the designer. 

IBBGG will be developed to allow complete 
generality in turtxxrachinery blade geometry 
specification . This provides a sasans to analyze and 
design all types of blades, and aft swept 

fan blades, axial and radial, exxp reisers, and axial 
and radial turbines. These features will be 
incorporated into the design/analysis system 
currently being developed at NASA Lewis. 

Either B-Spline data or NatiHJniform Rational 
B-Spline (NURBS) data will be used to generate 
surface coordinates for each blade row. This will 
allow the data to be ocopatible with CAD/ CAM 
systems. Instead of storing many blade coordinates, 
it is hoped that the above methodology will provide 
easier access to blade exordinates for manufacturing 
and as well as for finite element modeling for 
structural and dynamic analysis. Since IBBGG is 
currently in the development stage , formal 
discussion of the code and the specific methodology 
of the internal workings of the code will be 
provided in a subsequent report. 


In the middl e of figure 1, there is a block 
which serves as an interface between elements.^ This 
block ties together the grid, geometry, and flew 
solvers. As shown in figure 4, there are various 
(xepenents < x ■< men to all the major block elements. 
First, there is the geometry describing the 
turbccachinery problem cm is interested in solving. 
The geometry consists of blade and flcwpath 
information essential to the interactive grid 
generator, the visorus/ inviscid axisimmetrie 
design/analysis code, quasl-3-0 blade-to-blade 
calculations, and the full 3-D solutions. This 
geometry is also important to the interactive blade 
element package, because once the blades have been 
created and the data stored, the blades can easily 
be reconstructed for further analysis. 
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Figure 4. Interfacing Block for LERC Turbomochlnery 
Design/Analysis System 


Tne. grid and flow file information is 
essential to all flow solvers. The grid file can be 
generated by the interactive grid code or 
independently from existing grid codes. The flow 
file connects all the elements together. If an 
axisyircretric streamline curvature solution has been 
computed, this can be used to find the flow field 
from the VIADAC computation or as a starting guess 
for the 3-D flow solvers. The 3-D solution can be 
averaged axisynmetrically and used for the VIADAC 
computation. This is an important feature of this 
design/analysis system which allows closure between 
the full 3D and axisymretric approximation to the 
flow field. The grid file and flew file are written 
in a binary format carpatible with the IRIS and Cray 
Computers in a PLOT3D readable format, and can be 
submitted directly to each computer without pre- and 
post-processing these files. 

Another .important feature of the interface 
block is the "Blade Forces", The "Blade Forces" as 
described by Ariairezyk [19], are generally obtained 
from the 3-D solution, however, these forces can be 
obtained independently from quasi-3-D solutions or 
the forces can be modeled. These forces are 
important in order to ccnpute the work for a single 
or mil ti -blade- rew turbcmachinery problem with the 
design/analysis code. 

Finally, there are blade blockages, boundary 
layer blockages and blockage due to wakes or 
secondary flows. It should be pointed out that the 
boundary layers are cccputed if the axisymmetric 
flow code is run viscously. Wakes and secondary 
flew effects are to be modeled. These effects may 
also be included in the inviscid operation of the 
code in the form of a blockage. However, they are 
not essential to the operation of the code while 
running inviscidly. The blade blockage is always 
used while running the axisymretric confutation in 
the multi-blade-rcw or single-blade-rcw operation if 
a blade is present in the confutation. 

Viscous /Inviscid Axisymretric Desiqn/Analvsis Code 
(VIADAC) 

The third element of the design system is the 
flow code itself. VIADAC, the visoous/inviscid 
axisymretric design/analysis code, is formulated as 
an "average-passage" axisymmetric multi -blade- row 
flow solver. The code uses a discretized cell 


centered control volume approach with a feur-stage 
Runge-Kutta Scheme. Residual smoothing has been 
added which allows a larger Courant number and 
accelerate the convergence. Second and fourth 
difference artificial dissipation terms are used to 
stabilize the solution. The flow code incorporates 
a Baldwin-Lccax algebraic turbulence model for 
boundary layer confutations . Wall functions are 
automatical ly used when the grid stretching into the 
boundary layer falls outside the laminar sublayer. 
"Blade Forces" are used to simulate the effects of 
blade rows acting on the flew field. The spanwise 
inixing and secondary flew effects are currently 
being added to the flew solver. 

The flow solver can be operated independently 
of the system being developed at NASA lewis Research 
center. Although the current approach to the 
design of turtcrachines can perform the viscous flow 
confutation, the design system being developed will 
have the ability to use a streamline curvature 
calculation as a starting guess and run without the 
shear terms being confuted if desired. 

Currently the axisymretric part of VIADAC is 
operational , but the blade force and blockage terms 
have yet to be added. The modeling of the spanwise 
mixing and secondary flew effects are still being 
investigated. The axisymretric code has been tested 
on three different configurations and the results 
will be examined below. A publication describing 
the details of the system of equations and other 
modeling in the flow solver will be published in the 
future. 


and Discussion 


Several test cases were run to validate the 
VIADAC code. These cases are representative of 
typical axial turbcrrechirrery duct and axisyimretric 
wind tunnel body problems. Flow solutions were 
obtained from the NASA Rotor 67 Fan Test Program 
[201 , herein referred to as Rotor 67 Fan Test Bed. 
Rotor 67 Fan was tested at NASA Lewis Research 
Center as a single stage compressor with laser 
anemotretry data being acquired. There has been a 
great deal of valuable test data taken on Rotor 67 
Fan and are available to compare with to the 
solutions obtained by VIADAC. Eventually this 
involves including the force terms to simulate the 
multi-blade-row operation of the average-passage 
axisymretric code. 

The SR7 Spinner axisymretxic body was a 
configuration which was used during a propfan test 
in the wind tunnels at NASA Lewis Research Center. 
Experimental data were available can the axisymretric 
SR 7 spinner gearetry from an earlier test without 
the propfan blades attached to the spinner. The SR7 
spinner has seme interesting flew features along the 
spinner afterbody which test the axisymretric code 
capabilities along with the grid generation. 

A transonic bcattail body, tested by 
Shrewsbury [21], used to study the performance of 
air-breathing propulsive systems of supersonic 
aircraft operating at subsonic cruise conditions has 
provided an opportunity to study shcck-bcundary 
layer interaction with VIADAC. 


VXAEAC was operated without blade forces for 
all the cases to be shewn in this paper. First, the 
axisymretric grid for each case was generated using 
ITGGERC. In Figure 5, the grids for the SR7 spinner 
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NASA Bolor 67 Fan Test Bed Grid 162x61 Grid 



Figure 5. SR? , Rotor 67 Fan and Transonic Boattail Grid Geometries 


body, Potor 67 Fan and Transonic Boattail are shewn 
in blocks generated by TXSGSSC. Sfrtice that all the 
grids have different nuHibet of grid nod&s. The grid 
packing was varied from case to case to refine the 
boundary geometry where important flew features were 
likely to occur. The grid for the Rotor 67 Fan was 
blocked to incorporate the actual measure m ent planes 
upstream and downstream of the rotor. The grid for 
the SR7 Spinner body was packed around the spinner 
nose where large gradients in the flew field are 
expected and at the downstream end of the spinner 
where there is a step in the boundary. The grid for 
the boattail made! is packed on the afterbody, to 
allow accurate calculation of the transonic drag 
rise that occurs in this region. 

The computed surf ace static pressures and 
isentropic raach nunbers for the SR7 spinner body are 
shewn in comparison to the test data, in Figures 6 
and 7. These results were comforted assuming a free 
stream boundary condition along the upper grid line. 
All the results are presented in norriimensional 
form. The test data was acquired downstream of the 
spinner portion of the body. Static Pressure taps 
were located behind the spinner along the afterbody. 
The afterbody becomes cylindrical from z/IrO.1101 to 
z/Lf=0.4540 and ooriical from z/Xi=0.4540 to z/l>1.00. 
The surface geometry of the nose of the spinner and 
part of the nacelle were smoothed to ensure 
convergence of the flew solver. The actual spinner 
geometry has a very smal l step along the hub which 
was later incorporated into the grid. The solution 
obtained for that geemetry will be discussed later. 
The test data indicated a sharp drop with 'a sharp- 
rise in the static pressure just aft of the spinner. 



Axiat DisUfic* (2/L) 

Fl 8 uf. «• Pr.dlet.d 8ur<«e» St.tlc Pr...ur. lor 5*7 
Splnnar »©«*? 
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Figure 7. Predicted Surface Inaantropic Mach Numbers for 
SR7 Splrmar Bmty 


This is due to a sonic bubble along the surface of 
the afterbody as shown by the isentrcpic Mach number 
plot in figure 7. The ccrrputed results show a 
weaker shock at this point than that indicated by 
the data, however, the results closely ratch the 
data downstream and even where the geometry changes 
from a cylindrical to conic section, around z/L « 
0.4540, where a small rise in the static pressure is 
detected. The reason there is a weaker shock in the 
computation is due to smoothing of the boundary step 
ahead of the shock region. The sonic Nibble can 
clearly be seen along the hub surface in the static 
pressure and Mach number cxaitcurs shewn in figure 8. 

As mentioned earlier, the step for the SR7 
spinner was eventually added to the flcwpath and the 
solution cenputad. Particle traces shewn in figure 
9 along the wall near the step at z/L=0.0425 show a 
recirculation region and reattachment of the flow 
aft of the step. Another separated region existed 
upstream on the spinner body around z/L=>-0.030. 
This region was near the portion of the spinner 
where the blades would be attached. The small 
separation is evident without the blades on the 
spinner. Once the blades have been added to the 
spinner, the blockage # would increase, probably 
eliminating the separation. 

Rotor 67 Fan was analyzed because of the large 
experimental data base for both an isolated fan 
rotor and a full stage. Since no experimental data 
existed on the "duct only" operation, a second 
analysis was run on the duct using the same grid. 
For this purpose, PAKC2D [22] was used to perform a 
second set of carputations for comparison with 
results obtained from VIADAC. PARC2D is a canplete 
Navier-Stckes (2-D) flow solver in which the viscous 
tents can be selectively calculated with a thin- 


Static Pressure Distribution for SR7 Spinner Body 




Figure 8. Static Pressure and Mach Number Distributions 
For SR7 Spinner Body 
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Figure 9. Particle Traces Along the Endwall for the SR7 Spinner Body 


Mach Number Distribution From VIADAC 



Figure 10. Mach Number Distributions From VIADAC and PARC2D 
For Rotor 67 Fan Test Bed (Duct Only) 
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layer s inflation. The code uses a Beam and Warming 
approximate factorization algorithm. The algorithm 
is an implicit central difference scheme with a 
implicit artificial dissipation with a Jameson-style 
second-order term for improved shock capturing. A 
comparison of the Mach number contours for the "duct 
only” are shewn in figure 10, the results are in 
very good agreement. The surface static pressure 
distributions are shewn in figure 11. The results 


Predicted Surface Static Pressures 



Figure 11. Pr«dlet#d Surface Static Pressures From VIADAC 
and PARC2D for Rotor 67 Fan Duct 
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Figure 12a. Calculatad Radial Total Praaauras from VIADAC 

and PARC2D @ z/L » -0.061 for Rotor 67 Fan Duct 



Figure 12b. Boundary Layar Comparison of Total Pressure from 
VIADAC and PARC2D @ z/L - -0.061 for Rotor 67 




from the two codes are identical. Total pressure 
profiles at z/I>H3.061 just upstream of the leading 
edge of the rotor shown in figure 12a are in fairly 
good agreement. As shown in Figure 12b, PARC2D 
produces a very snail overshoot in the total 
pressure near the edge of the bcundary layer while 
no overshoot in the distribution is detected frem 
the VIADAC computation . The artificial dissipation 
is turned off with the present code by locating the 
edge of the boundary layer in the oorputation, while 
PARCS D distributes the artificial dissipation to the 
walls. VIADAC used approximately 300 epu sec on 
the Cray YMP ard converged in about 3000 cycles. To 
completely eliminate the total pressure overshoot in 
PARC2D, PARC2D would require about 10,000 or more 
cycles to ensure proper convergence had been 
achieved. Eventually, this case will be run with 
VIADAC including the blade forces for the rotor 
only, then the stage. 


Finally, solutions were obtained for an 
axisyrnoetric boattail geometry model [21] , which had 
a boattail angle of 15 degrees and F/EKL.0, at three 
different Mach numbers to investigate the transonic 
drag rise and shock formation on the afterbody. The 
confuted pressure coefficients for the three Mach 
lumbers , Mach - 0.56, 0.70, and 0.90, are compared 
with wind tunnel data in figure 13 ard are in very 
good agreement. The flew accelerates ante the 
af terbo dy at z/If0.9 ard forms a weak shock at a 
free stream Mach number of 0.90. This is more 
evident in the Mach contours shewn in figure 14. The 
weak shock forms around z/D=0.94 and perhaps a very 
small separation region forms near the exhaust 
region around z/Ip 1.00. It is not clear from the 
test data as to whether this separation is large, 
small, or nonexistent, hewerver, the code was able to 
predict the shock location fairly well. 
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Figure 13. Calculated Surface Preasurea from VIADAC for a 
Transonic Boattail @ 0.56,0.70 and 0.90 


A design/analysis system for 

viscous/ inviscid flews through turb^ohinery ^ 
currently being added to the flew codes at NASA 
Lewis Research Center. Grid generation is always a 
crucial part in any analysis using a Navier-Stotes 
flow code. TXGGERC is able to grid sophisticated 
aeanetry shapes and versatile enough to create the 
STdLS packing nec^ry for carpitrng 
implex flows. This is a crucial step 
skmificartt capability in solving the flows in 
turbcnachinery. As dencrstoted^e^ilihy 
grid real geometries, such i i as the step ^ 

available in the past. 

Solutions obtained from the axisyntnetric 
^ solver/ VIADAC, have shown to be in 

obtained from PAPC2D. The SK7 spinner and 

Sc boattail exiles established the <£uity 
of the codes to solve and analyze complex 
aLeSries The ability of the code to predict 
separation and shock boundary layer interaction^ 
critical, sinoe the code is to be run in all types 
of regimes. 



Figure 14 . Computed Mach Number Contours 
for the Transonic Boattail AfterBody 
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The ccinp arisons between VTADAC and PAKC2D were 
shown to be in good agreement, The solutions 
cbtained for Rotor 67 fan are only a starting point 
in the development of VIADAC. currently, blade 
forces for the rotor are being added to the 
solution. Future efforts will include adding a 
more sophisticated model for spanwise mixing and 
turbulent diffusion for each blade rw. 

The grid generator is cccplete and has already 
proved to be a valuable tool in investigating real 
physics. The flow solver has provided excellent 
agreement to "duct only" flew fields and is 
currently being upgraded to include the forces 
necessary to compute the flews in railti-blade-rcw 
turbcmachinery problems. Future efforts will 
include the completion of the surface geometry 
generator for the design of the blades. 
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